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We present a new method to detect Fermi surface instabilities for interacting systems at finite 
temperature. We first apply it to a list of cases studied previously, recovering already known results 
in a very economic way, and obtaining most of the information on the phase diagram analytically. As 
an example, in the continuum limit we obtain the critical temperature as an implicit function of the 
magnetic field and the chemical potential Tc(/i, h). By applying the method to a model proposed to 
describe reentrant behavior in SrzRu207^ we reproduce the phase diagram obtained experimentally 
and show the presence of a non-Fermi Liquid region at temperatures above the nematic phase. 

PACS numbers: 75.10 Jm, 75.10 Pq, 75.60 Ej 



I. INTRODUCTION 

Landau's theory of the Fermi liquid (FL)^ ^ is the 
standard description of weakly interacting fermionic 
systems. The basic assumption is that, at the micro- 
scopic level, there is a statistical ensemble of elementary 
excitations that behave as fermionic quasiparticles with 
well defined momenta and energies. These quasipar- 
ticles are assumed to have long lifetimes, and their 
distribution function shows a step at finite momentum 
defining the Fermi surface (FS) in momentum space. 
In this picture, low energy excitations are described 
by infinitesimal deformations of the FS. When one of 
such distortions leads to a reduction in the free energy, 
we have a breakdown of the FL description due to 
what is known as a Pomeranchuk instability^. This 
kind of instabilities has been intensively studied during 
the last fifteen years. In particular, broken symmetry 
phases arise due to electron-electron interactions, within 
the Landau-Fermi liquid description^ In some 

cases, such instabilities lead to nematic phaseJ^^M^. A 
comprehensive review about nematic Fermi fluids with 
an extensive list of references could be found in Ref. [31] . 

In previous works we developed and applied a general 
method to detect Pomeranchuk instabilities. With it, 
we were able to obtain the instability regions for several 
systems. While the first work^^ was mainly advocated 
to study cases involving simple interactions (as 5-wave 
and d-wave interactions) and zero-temperature square 
lattice systems, the method was extended to incorpo- 
rate temperature dependence (mainly within a low tem- 
peratur e exp ansion approximation) and flavor degrees of 
freedonP^^. Also, application to systems with differ- 
ent lattices, as the honeycomb case has been successfully 
accomplished^^. The approach consisted in decomposing 
the deformations of the FS in a set of basis functions or 
"channels", and then studying the contribution of each 
such channel to the free energy. 

It is the aim of this work to present a more efficient 



method to diagnose FS instabilities. With this new 
approach, by making use of the symmetries of the 
interaction function, we are able to find the FS instabil- 
ities in a much more economic way, as compared with 
previous approaches. We can now obtain the unstable 
regions avoiding much of the time-consuming numerical 
calculations that were faced previously. Furthermore, 
with this method there is no need of any extrapolation 
to higher channels to determine the phase boundaries, 
which in many cases can be obtained analytically. In 
addition, the present method allows to study finite 
temperature problems without the need of a low tem- 
perature expansion. The method is discussed in detail 
in Section [Til 



In Section III we apply this method to several cases for 
both zero and finite temperature. There we show that 
applying the present method to the previously studied 
cases, the reported results are recovered. Moreover, for 
much of the cases discussed in this section, we get the 
instability condition as an analytic expression of the 
parameters involved. For instance, when we study the 
stability of a finite temperature system in an external 
magnetic field with a dispersion relation corresponding 
to the continuum limit, we are able to find an implicit 
relation for the critical temperature as a function of the 
applied magnetic field and the chemical potential. 



In Section IIVI we consider a d-wave interaction func- 
tion. This kind of interaction has been proposed 
to describe nematic Fermi liquid^^ and shown to be 
relevant in the discussion of the reentrant behavior 
manifested i n ultr a- clean samples of Sr2,Ru20j at low 
temperature j^^ l ^A There, as the external magnetic field 
is increased, the compound goes from an isotropic Fermi 
liquid phase to a nematic one, and then back to the 
isotropic phase. The intermediate phase forms a dome- 
shaped region in the T vs. h phase diagram. A review 
of the physical properties of Sr2,Ru20j can be found in 

By studying this model within a mean field approxima- 
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FIG. 1. (Color online) Schematic phase diagram, showing 
both Fermi liquid (FL) and no n- Fermi liquid phases. The con- 
tinuous (dashed) lines correspond to continuous (discontinu- 
ous) phase transitions. Gray dotted lines correspond to con- 
tinuous transitions between isotropic and nematic FL phases 
that arise within our stability analysis but which are shown 
to be preempted by a discontinuous transition. See Section 
HVlfor details. 

tion, Yamase and co-worker^^ obtained - for some values 
of the parameter's space - a wide region of stability for 
the isotropic Fermi liquid and the characteristic dome- 
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Here e^{k) is the dispersion relation that controls the 
free dynamics of the system, and the interaction function 
/"^(k, k') can be related to the low energy limit of the 
two particle vertex. We will decompose the interaction 
function as 

N 

r''(k, k') = f/'^d?(k)4(k') (2.2) 

where {d^}i^i represents a given basis of the space of 
functions on momentum space. For most of the inter- 
actions used in the literature, such basis can always be 
chosen in such a way that the sum has a finite number 
of terms N . 

With this at hand, we will prove that the variation of 
the Landau free energy can be rewritten as 

N 

Sn = ^c'c^Mij (2.3) 



shaped region for the nematic Fermi liquid. They also 
obtained the observed behavior for the phase transitions 
between both phases (continuous transitions for the up- 
per temperature part of the dome and discontinuous for 
the lower temperature part). Since then, it has been well 
establishecP^ that this approach is a good starting point 
in the construction of a phenomenological description of 
the observed behavior in Sr^Ru207. 

A major problem of the mean-field approach, when 
used alone, is that it fails to describe the non Fermi 
liquid region of the phase diagram: the model seems 
to indicate that the isotropic Fermi liquid phase is also 
present there, even for higher temperatures, above the 
dome-shaped region. Our goal is to improve the mean- 
field approach, by combining it with the Pomeranchuk 
stability analysis developed in the present paper, to 
obtain a phase diagram schematically shown in Fig.l. 

Section W\ contains the conclusions. 



II. THE METHOD 

According to Landau's theory of the Fermi liquid, the 
change in the Landau free-energy = E — TS — /J^^Na 
can be written as a functional of the change Jn"(k) in 
the equilibrium distribution function at finite chemical 
potential /i" for each particle species a e [1 . . . Nf]. To 
first order in the interaction. It reads 



(2.1) 



I 

where {c^}ie[i,..,N] are arbitrary real numbers param- 
eterizing the deformations of the occupation numbers 
(5n"(k), and we have defined the "stability matrix" Mij 
in the form 

N 

M,j = (o^,|d,) + ^/7^^(4|^z)(^id^,) (2.4) 

kl 

in terms of the positive-definite bilinear forms 

m)= f rf'fc^(e"(k)-Ai")F"[^"-e"(k)]r(k)r(k) 

im = f d2fc^F'[^"-e"(k)]r(k)r(k). (2.5) 

The key idea in Pomeranchuk's analysis is to characterize 
deformations (5n"(k) that would lead to Sft < 0, then 
pointing to an instability of the system. Then, according 
to ( |2.3[ ), in order to have a positive energy for all sets of 
{c^}ie[i,..,N]j the matrix Mij has to be positive-definite. 



Sn = J2 yA(£"(k)-M")<5n„(k) + i^^ Jd\Jd\'r^{k,k')Sn^{k)Snp{k') 



a 13 
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A simple way to check positivity is to verify that ah the 
minors of the matrix Mij are positive, i.e. 



Min 



I Minors (M) I 



> 



(2.6) 



ie[i,N] 



It has to be pointed out that the calculation of the 
stability matrix involves only a finite number 2N'^ of in- 
tegrals, which in most cases are straightforward to com- 
pute, and in some of the problems studied in the litera- 
ture can be even found analytically. 



This is our central result: to find the Pomeranchuk insta- 
bilities of a given system, one has to compute the corre- 
sponding stability matrix (2.4), evaluate its minors which 



are functions of /i", T and Uij^ and find the region in 
phase space in which they become negative. 



Proof of the result 

In order to prove the result presented above, we insert 
the decomposition of the interaction function (2.2) into 



the variation of the Landau free energy (2.1), to get 



'^^ = E/ d\{e"{k)-ii'') 5na{k) + / d\df{k)Sna{k)^ ( J d\'d^ {k')Snp{k') 



(2.7) 



The next step is to parameterize any given deformation of 
the occupation numbers with the help of a new function 
(5"(k), as 

(5n"(k) = F[/x"-£"(k) +,5"(k)] -F[/u"-e"(k)] = 

~ F' [/Li"- (k)] 5" (k) + ^F" [/Li"- (k)] .5" (k) ^ 

(2.8) 

where F[x] is the Fermi distribution 
Fix] = ^ 



(2.9) 



and a prime C) denotes a derivative with respect to x. 

With the help of eqs. ( |2.8[ ), we can rewrite the varia- 
tion in 50. to lowest order in (5"(k) as 

5^=\ [m + Y.U'^{d^\5){d,\5)\ (2.10) 



The proof is as follows: given the basis {d^}i^x we sepa- 
rate it into two parts {d^}ie[i,...,N] U {(if }ie(/_[i,...,Ar]) 
and replace the second factor by another set of ba- 
sis functions {df }i^(/_[i^...^7v]) to which we have sub- 
tracted the part non-orthogonal to {df n]i namely 
Jf (k) = df (k) - Y.i^^^,...,N]{di\di){d,\d,)-^{]^). By a 
similar argument, we can further decompose 



N 



<5J(k) = ^cX"(k)+5!^(k) (2.13) 



where (5"_|_(k) is defined as satisfying 



Vi: (<5ii|di) = 



(2.14) 



where the bilinear forms were defined in (2.5). If we now 



decompose 8'^ in the same basis {d^}i^i as the interac- 
tion function, namely 



N 



(5" (k) = c'df (k) + (5J (k) (2.11) 



Recalling the previous condition in 6_\_ (k) we get 



N 



Yc'{d,\dj)^{S^^\dj)=0 (2.15) 



it is easy to prove that we can always choose (5J(k) as 
satisfying 



V, : {S^\di)=0 



(2.12) 



which implies that the c* are not independent parameters, 
but they are given in terms of the c*'s. Replacing into 
the expression (2.10) for energy we obtain 



AT / / AT 

= 9 E (^^1^^-) + E u''{dk\d,){di\dj) 



kl 



+ 2c'c^{di\dj)+c'(J{di\dj) J + - {5^^\5^^) (2.16) 
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By comparing the above expression with equation (2.10), 



we see that the last term represents the fluctuation of a 
free fermion with Uij = when perturbed with ^xx(k). 
This implies that this last term is always positive and 
does not contribute with instabilities to the above sum. 
Then we can safely put ^±±(k) = 0. This in turn implies 



that, according to eq.(2.15) and if the matrix {di\dj) is 
invert ible, then the coefficients P have to vanish. We 
then obtain 



iv / N 

5Q = -Y^c^c^ {d,\dj) + Yl U^\dk\d,){di\d^ 



kl 



(2.17) 



In conclusion, the energy can be written as the quadratic 
form (2.3), and we have proved our main result. 



Notice that the factor {di\dk) correspond to the stabil- 
ity matrix of a zero temperature Fermi system without 
interactions, i.e. a Fermi gas. Since such Fermi gas is 
obviously stable, then by the considerations of the pre- 
vious section, its matrix {di\dk) is positive definite. The 
condition of positivity of Mij is reduced to the positiv- 

ity of the "zero temperature stability matrix" . In 
other words, at zero temperature the Fermi liquid will be 
stable as long as 



Min|Minori(M(^))| 



> (3.4) 

' *e[l,A^] 

In what follows we apply this condition to some simple 
cases 



Spinless case 



III. APPLICATION TO KNOWN EXAMPLES 

In this section we will apply the method presented 
above to systems whose Pomeranchuk instabilities have 
been pr eviously s tudied in the literature by alternative 
method^^ ^ ' ^^ ' ^^'^^*^ . We will start by applying it to the 
zero temperature spinless case, and progress later to more 
complicated situations. 



III.l. Zero temperature 

In the zero temperature limit, both bilinear forms 
(2.5) contains a Dirac ^-function in the /i"— e(k) vari- 
able. Indeed, when T ^ we have F'[x] 5{x) and 
—xF''[x] -x3'{x) = S{x). This implies that the bilin- 



ear forms coincide, and depend only on the values of the 
integrand evaluated at the Fermi surface 

(3.1) 

This integral can be further simplified by parameterizing 
the Fermi surface of species a with the help of a param- 
eter 5" G 5^ as k = k(s"). We gelP 

(3.2) 

An immediate consequence of this is that all the c* exci- 
tations vanish. The stability matrix reads 



N 



N 



N 



kj 



(3.3) 



At zero temperature and in the absence of spin, the 
index a is not present and the sums in the bilinear forms 
(3.1) collapse to a single term. The resulting expression 



is particularly simple 



^*ML«k(,)Wk(,)) (3^.) 

We will start by studying the simplest possible case of 
an interaction containing a single function d{k.). In other 
words, we will chose an interaction function that can be 
written as 



/(k,k') = /7d(k)d(k') 



(3.6) 



In this case the zero temperature stability matrix con- 
tains a single entry, and the stability condition (|3.4[) reads 



U{d\d) > 



(3.7) 



with {d\d) as in (3.5). 

As the first concrete example of the above simplifica- 
tions, we can use the method to study the Fermi liquid 
in the square lattice that was previously investigated by 
a different method in [32]. In a square lattice, the Fermi 
surface is defined by 



£:(k) = —2t (cos kx + cos ky) = ji 



(3.E 



where t is the hopping parameter, and ji the chemical 
potential. It can be parameterized as 



kx = arccos I — ^ — \- s \ 



ky = arccos (^———s 



At 



M 

4t 



< 5 < 1 



where —1 

bilinear form (|3.5|) reads 



At ' 



(3.9) 

With this at hand, the 



(3.10) 
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FIG. 2. Instability region for a Fermi liquid with an s-wave (a) 
and a d-wave (b) interaction in the square lattice. Note that 
both plots coincide with the extrapolation to higher channels 
of the plots presented in Ref . [32] . 



Now in the case of an s-wave interaction in the square 
lattice, we have 



/(k,k') = C/, 



(3.11) 



that can be decomposed in the form we need by using a 
single function (i(k) given by 



d(k) 



(3.12) 



In this case the integral ( |3.1Q ) can be performed explic- 
itly, an we get an analytic result for the stability condi- 
tion in parameter space 



16 [/ 



4t ■ 



Re 



K 



(4^ - m)^ 
(4t + m)^ 



> 



(3.13) 



where K{- ■ ■) denotes the complete elliptic integral of the 
first kind. With this expression we immediately get the 
phase diagram of Figj2|a). Note that the plot coincides 
remarkably with the extrapolation to higher channels of 
the plot presented in Ref. [32 . 

If instead we have a d-wave interaction, namely 

/(k, k') = —U (cos kx — cos ky){cos k'^ — cos k'y) (3.14) 

that can again be represented in terms of a single basis 
function (i(k), given in this case by 

d{k) = coskx — cos ky = 2s (3.15) 

after performing the integral (3.10) to obtain {d\d) we get 
for the stability condition M^^^ > 



4[/ 

t2 



(4t + /i)Re 



K 



(4t-M)^ 
(4t + M)^ 



(4t - fi)' 
(4t + /i)^ 



> 
(3.16) 

where E{- ■ ■) represents the complete elliptic integral of 
the second kind. The resulting phase diagram is shown 
in Fig 2|b). Again the diagram coincides remarkably 
with the extrapolation to higher channels of the plots 
presented in Ref. [32] . 



Now we can progress into more complicated situations. 
For example we can consider non-square lattices, and/or 




V 

0.75 
0.5 
0.25 


-0.25 
-0.5 



(&) 



1.5 2 2.5 




FIG. 3. Instability regions (shaded) for a Fermi liquid with 
the interactions given by (3.17) in the honeycomb lattice. 
Planes ^; = (in fig. (a)) and /i = 0.95 (fig. (b)). Both 
plots seem to recover the extrapolation to higher channels of 
what was reported Ref. [35] . 



interactions with multiple components, in which case the 
zero temperature stability matrix would contain multiple 
entries. 

As an example, in Ref. [35] a honeycomb lattice was 
studied. There, on site interactions with strength u and 
nearest neighbor interactions with strength v were con- 
sidered. In its factorized form, the interaction function 
reads 

7 

/(k,kO = ^/7M,(k)d,(kO, (3.17) 

i=l 

where the coefficients are = {2u-\-3v)/{4:7t)'^ and W = 
—v/{4:7t)'^ for i = 2 to 7, and the functions (ii(k) are 



rfi(k) 
d2(k) 

d6(k) 



(3.18) 



cos (k.(5i 
cos (k.Ss 
sin (k.(52 - 

V3\ 



-<P) 



d3(k) 
d5(k) 
rf7(k) : 



: COS (k.(52 
: sin (k. - 
sin (k.c^s - 



with = a{\, ^), 32 = a(^, ^) and S3 = a(-l,0). 

Introducing this information into the expression of our 
zero temperature stability matrix, it becomes a 7 x 7 
matrix given by 



io) 



Sij ^U'{di\dj) > 0. 



(3.19) 



Using it in (3.4) we get the phase diagram of Figjs] 
Again, notice the coincidence with the results obtained 
by different methods. 



Spinful case in an external magnetic field 

A further step would be the inclusion of spin, which 
allows for the coupling to an external magnetic field. We 
will limit our study to the case of a simple interaction 
that can be written in terms of a single (i"(k) as 



r^(k,kO = /7d"(k)d^(kO 



(3.20) 
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FIG. 4. (Color online) (a)-(d) Instability regions (shaded) for 
external magnetic field going from /i/t = to /i/t = 3. The 
plots seem to recover the extrapolation to higher channels of 
what was reported Ref.[34 . 



In the case of a spin antisymmetric interaction in the 
square lattice with nearest neighbor hopping, the sim- 
plest example of the kind studied in [34] is given by 



d"(k) 



(3.21) 



with a = ±1. With this interaction and recalling the 
results of the previous section for the square lattice, we 
immediately get from the stability condition M^^^ > 
the requirement 



E 



16/7 



4t-\-/j.-\-ah 



Re 



^f !r""i!n i>o (3-22) 



where 2h is the applied magnetic field. 

A somewhat more complicated interaction studied in 
[34] is the one given by 



(3.23) 



Remarkably, the zero temperature stability matrix for 
this interaction is given by the same formula ( |3.22 ), with 
U replaced by /7/2. The result coincides with the spin- 
split instability found in [34 , and in particular, as it can 
be seen in Figj4] reproduces the re-entrant behavior re- 
ported there, that was observed in experiments. 



III. 2. Finite temperature 

With the training of the previous sections, we can now 
proceed to study the finite temperature case. Now, the 
bilinear forms introduced in (2.5) are no longer equal. 



(b) 



V 
-0.1 





FIG. 5. (Color online) Instability regions (shaded) for s-wave 
interaction function, (a) 3D phase diagram showing unstable 
regions in the space of parameters /i, JJ and T. (b) Instability 
regions for kT — 0.01,0.5 and 1. (c) Instability regions for 
-0.2. 



Spinless case 

With a simple interaction in terms of a single (i(k) as 
in (3.6), we get the stability matrix as 



M = {d\d)^U{d\df >0 



(3.24) 



In the case of an 5-wave interaction (3.11 )-( 3.12 ), the 



integrals in the {d\d) and {d\d) brackets have to be eval- 
uated numerically. We get for the phase diagram the one 
displayed in FigjS] 

Again, note that the stability region coincides with 
the previously published results obtained by a different 
method in Ref.[33]. In particular, in that work an ex- 
pansion in powers of kT/fi was needed, while the present 
calculation is exact to any order in this parameter. 



Spinful case in an external magnetic field 

To test our method on a spinful system at finite tem- 
perature, we will use the simple dispersion relation cor- 
responding to the continuum limit, i.e. 



and the stability matrix may be more involved. 



1 

2m 



+ ky 



(3.25) 
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FIG. 6. (Color online) Instability regions (shaded) for differ- 
ent fillings (from light to dark coloured lines) /i = 0.1, 0.2, 0.5 
and interaction strength 2nmU = 0.1 in fig(a) and 2nmU = 
0.15 in fig(b). The boundaries of such regions indicate the 
critical temperature as a function of the external magnetic 
field. 



in other words, no underlying lattice is considered here. 

In the case of the spin-antisymmetric interaction given 
by eqs. ( 3.20 )-( 3.21 ), the integrals in the {d\d) and {d\d) 



brackets are easily evaluated by turning into polar coor- 
dinates in the ky plane. We get 

{d\d) = 27rm^(F[/i + a/i] - {jj. ^ ah) F' [jj. ^ ah]) 

a 

{d\d) = 27rm ^ F[fi + ah] (3.26) 



from which the stability condition M > can be easily 
obtained 



27rmU 



^(F[/i+a/i]-(/i+a/i)F'[/i+a/i]) 



(3.27) 

With this expression, we can plot, e.g.^ the critical mag- 
netic field he as a function of /i,T for any given [/, as 
shown in Figj6] 



IV. APPLICATION TO Sr3Ru207 



In this section we apply the present method to the 
model HamiltoniarW 

H = ^l (e<'(k)-M-a/i)na(k) + ^^r''(k,k')n„(k)n;3(k) | 

(4.1) 

where the dispersion relation is that of a 2-dimensional 
square lattice with first and second neighbors hopping. 

e^Ck) = —2 (t(cos kx + cos ky) + 2t' cos k^ cos ky) (4.2) 

The interaction function takes into account only forward 
scattering and has the form 

/"^(k,kO = -/7e"^ {coskx- cos ky) {cos k'^ - cos k'y) , 

(4.3) 



where e^^ is a matrix with all its entries equal to 1. As 
discussed in section |n| this interaction is separable with 
the function (i(k) taking the form of a d-wave form factor 

(i"(k) = e"(cos/c^ - cos ky) (4.4) 

where is a vector with all its entries equal to 1. 

IV. 1. Mean field calculation 

Within the mean field approximation, the free energy 
reads^^: 



F^Yl (^'(^) -f^-o^h- £«(k)) n«(k) 
^ ^ r^(k,k0nk.nk^;3- 



k.Q; 



27V 

fcBT^ln(l + e-^"('^)/'=-^). (4.5) 



By minimizing it, we find the self-consitency equation 
s«(k) = e%k) -^-ah- ^c^"(k)7?. (4.6) 



where use has been made of eqs. (4.3) and (4.4). In this 
expression we have defined 



k,/3 ^ 



(4.7) 



For each set of parameters {/i, T, /i, /7, t'}^ 77 is a real 
number to be found self consistently. Indeed, by writing 
explicitly the dispersion relatio n g^ (k) and the d-wave 
form factor (i"(k) in expression (4.7), we get 



Sa{k) = - (2t + ^) COskx - (2t - ^) 



COS k., 



-\-4:t' cos kx cos ky 



ah, 



(4.8) 



Here we see that whenever 77 vanishes, the system has 
7r/4 rotational invariance. On the other hand whenever 
?7 is different from zero, isotropy is broken and we have 
a nematic phase^^. For that reason, we will call rj our 
"order parameter" from now on. 

First of all, notice that the solution r] = (the isotropic 
Fermi liquid solution) is consistent at any point of the 
parameter space. 

In addition, o ne a ls o fin ds numerically the existence of 
solutions of eqs. (4. 6)- (4. 7) with r] ^ 0. We have explored 



the h vs. T plane for t = 1^ f = 0.35t, U = /i = 1, 
and the results obtained for the order parameter 77 as a 
function of temperature and magnetic field are shown in 
Fig[7|^a). As expected, these results coincide with those 
presented in Ref. [39] , implying the presence of a nematic 
Fermi liquid phase inside a bounded region of the h vs. 
T plane. This region is shown in FigjT^b). 
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Both unstable 





FIG. 7. (Color online) (a) Order parameter as a function of 
magnetic field h and temperature T. (b) Result obtained for 
U — 1 and 11 — 1. The solution with 77 = exist at any 
point of the parameter's space. On the other hand, the ?7 / 
solution only exists inside the grey dome-shaped region, 
(a) (h) 





h 



h 



FIG. 8. (Color online) (a) According to the mean field anal- 
ysis the isotropic phase [i.e. r/ = 0) is defined in all the T vs. 
h plane. The instability check performed with the extended 
Pomeranchuk method, shows that it is unstable in the blue 
region and stable in the rest of the plane, (b) The nematic 
phase (i.e. 77 7^ 0) is defined only inside the dome shaped 
region. According to the extended Pomeranchuk method, it 
becomes unstable in the green region. 



IV. 2. Extended Pomeranchuk analysis 

With the numerical values obtained for the order pa- 
rameter 77 inserted into the dispersion relation, we can 
use our extended Pomeranchuk method to check the sta- 
bility of the isotropic and nematic phases at any given 
point of the parameter space at which they are defined. 

The isotropic phase is a solution of the self consistency 
equations at any point of the parameter space, and us- 
ing our method we find that it is unstable in the region 
depicted in Fig. (sja). On the other hand, the nematic 
phase is defined in the grey region of Fig. [Tfb) and we 
find that it is unstable at the contour of the dome as is 
indicated in Fig. [sj^b). 



We can now put together these results to construct the 
phase diagram shown in Figj9j which is the main result of 
this section. There we can distinguish the two different 
phases: the isotropic Fermi liquid phase, corresponding 



Isotropic 




Isotropic 



h 

FIG. 9. (Color online.) Phase diagram: In the white and 
green regions, only the isotropic Fermi liquid phase is stable 
while in the red region the only stable phase is the nematic 
Fermi liquid. Inside the grey regions both phases are stable 
and finally, in the blue one, neither of both Fermi liquids are 
stable. 



to white and green regions and the nematic Fermi liquid 
phase in the red region. Remarkably, by increasing the 
temperature it is possible to reach the blue region, where 
both (isotropic and nematic) Fermi liquids are unstable, 
indicating the breakdown of the Fermi liquid description. 
In this way, by complementing the mean-field analysis 
with the extended Pomeranchuk analysis, the non-Fermi 
liquid region of the phase diagram is obtained. 

Yet another interesting feature of the obtained phase 
diagram is the existence of regions where both phases 
are predicted to be stable. These regions are indicated in 
gray m Fig. [9| As our method could only predict instabil- 
ities associated with continuous phase transitions, wher- 
ever we detect a possible coexistence of different phases, 
we must compare their free energies in order to identify 
the actual phase, that would be reached via a discontin- 
uous transition, as discussed in Ref. |39l 



IV. 3. Discontinuous phase transitions 

In the presented phase diagram, we still have to de- 
termine which of the phases survive in the coexistence 
regions. To do that, we will use the mean field expres- 
sion for the free energy, eq. (4.5), replacing the ansatz 



(4.6) for £a(k) and using the d-wave interaction expres- 
sion (|43|. We get 



2N 



-/3£o(k)l 



(4.9) 
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F (c) 



FIG. 10. (Color online) (a) Zoom of the phase diagram 
showed in Fig|9] near one of the coexistence regions, (b) Free 
energy as a function of the order parameter 77, at T=0.08. 
There we see a continuous phase transition. When the tem- 
perature is lowered we have a discontinuous phase transition, 
as is shown in (c). 



We can now expand this free energy as a power series in 
the order parameter, obtaining 



F(7?) ^ F(0) + la2n 



4! 



04 ?7 



(4.10) 



Now with the help of the coefficients a2, a4, etc. we can 
plot the free energy in the coexistence region. The results 



are shown in Fig 10 We see that at low temperatures we 
have a discontinuous phase transition while, when tem- 
perature is increased, the transition becomes continuous. 
We also note that in the coexistence region the minimum 
of the free energy is achieved for 77 7^ 0, implying that 
the nematic phase is the one that prevails. 

In this section we have shown that the extended 
Pomeranchuk method is a useful tool to complement the 
mean-field description of the low temperature phase di- 



agram of Sr2,Ru20j^ reproducing both the isotropic and 
the nematic FL phases, as well as the non Fermi-liquid 
phase. 



V. CONCLUSIONS 

To summarize, in this paper we have presented a more 
efficient method to diagnose FS instabilities. Within this 
approach, by parameterizing appropriately the Landau 
interaction functions, FS instabilities are detected in a 
much more economic way as compared with previous ap- 
proaches. We can now obtain the unstable regions ana- 
lytically in most of the cases, and in those cases in which 
numerical evaluation is unavoidable, the numerical ef- 
fort needed is minimal. In addition no extrapolation to 
higher channels is needed to determine the phase bound- 
aries. More importantly, the present method allows to 
study finite temperature problems without resorting to a 
low temperature expansion. 

We have applied the method to previously studied 
cases and recovered all known results, while instability 
conditions are obtained as analytic expressions of the 
parameters involved. In particular, for a continuum 
model at finite temperature in an external magnetic 
field we have found an implicit relation for the critical 
temperature as a function of the applied field and the 
chemical potential. 

Finally, we have applied our method to a (i-wave in- 
teraction function in the square lattice, w hich has b een 
proposed to describe nematic Fermi liquids^ *^^ * ^^*'^-'^ *. By 
combining a mean-field analysis with the Pomeranchuk 
stability analysis developed in the present paper we re- 
produced many features of the experimentally observed 
phase diagram^^. 
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